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The relaxational modes and aging of the Glauber dynamics of the mean-field spin-glass 
model, SK model, are studied by a numerical diagonalization technique and Monte Carlo 
simulations. We found that the aging process of the model is understood as hierarchical 
growth of quasi-equilibrium domain in the phase-space. 

§1. Introduction 

The fascinating features of aging effects of spin-glasses revealed by experiments!^ 
have arose niuch theoretical interest in recent years, which includes various phe- 
nomenologiescPtlP, analytic predictionstS and numerical simulations. EP EP In order 
to seek for concrete information to provide theoretical base, we study the relax- 
ational modes of the Glauber dynamics of the mean-field spin-glass model, namely 
Sherrington-Kirkpatrick (SK) modelEP of finite sizes at temperatures below T c (spin- 
glass transition temperature). By numerically diagonalizing the transition matrix of 
small system sizes, we obtain a spectrum of relaxational modes and analyze their 
properties to get insight into the mechanism of the aging process. We also discuss 
the data of Monte Carlo simulations we have performed on larger system sizes. 

§2. Model 

We study the SK model, whose Hamiltonian is given as, 

H = s ^J i jO- i o-j - h^2,<Ji. (2-1) 

ij i 

where cij's are Ising spins, {Jij} are independent random Gaussian variables with 
zero mean and the variance (N — 1) _1 and h is the external magnetic field. 

We denote the probability that a spin configuration {erf} is realized at time t 
as p a (t). The time evolution of the probability p(t) = (pi(t),p2(t), . . . ,p 2 N(t)), is 
detemined by the master equation, 



d 
~dt 



p(i) = -r P (t) with -Y aP = W aP - 8 aP E 7 W 7/3 . (2-2) 



The transition probability Wq,^ is chosen as the heat-bath type, W Q/ 3 = ^(1 
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af tanh(T _1 (^j ■ Jij°~j + h))). The unit of time corresponds with 1 Monte Carlo 
step/spin (MCS) used in the standard Monte Carlo simulations. 

As usually done, it is convenient to introduce a matrix r a p defined through the 
similarity transformation T a p = Pcq^f a pP\l^ . The matrix f a p is a real-symmetric 

matrix and can be diagonalized as T a p = 2~2i=i u a{l)z{l)up{l) where u a (l) is the a-th 
element of the ortho-normal eigen- vector corresponding to the eigen-valuse z(l). We 
label the eigen-values in the order = z(l) < z{2) < ...z(2 N ). The first mode 
corrsponds with the equilibrium distribution and the eigen-vector is known exactly 
as u a (l) = \Jp°a. In our analysis, we perform numerical diagnalization (Householder 
method) to obtain the whole eigen- modes of system sizes N = 8, . . . , 12. 

Once the eigen-modes are obtained, G a p(t), the transition probability to go from 
(5 to a in time t, can be wirtten as 

2 jv 

G a p(t)=p^r a (l)exp(-z(l)t)rp(l) with r a (l) = u a (l)/y/j&. (2-3) 

i=i 

§3. Analysis 

3.1. Aging after Temperature Quench 

Now we consider aging process after rapid temeprature quench from the infinitely 
high temperature down to a temperature below T c . For this purpose, we choose the 
initial condition for the master equation as p a (0) = 1/2 N , which means we choose 
the initial spin configuration at random. Then the system ages as 

P 

by the propagator G a p(t w ) of a certain temperature T < T c . 

In order to measure the extent to which the system approachs the equilibrium 
with a given i w , it is convenient to introduce the following indicator of aging, 

Ra(t w ) = Pa(t w )/ P c a q = £ r a (l) exp(-z(/)t w ) ]T rp(l)pp(0). (3-2) 

l (5 

In Fig. 1 we plot an example of the evolution of the indicators at spin con- 
figurations of different energy minima. Suppose that the indicators of a pair of 
spin-configurations, say a and (3, satisfy the equality with certain characteristic time 

t* p] 

R a (t w ) = Rp(t w ) fortwX^. (3-3) 

Then we say that the two configurations are in quasi- equilibrium, because the relative 
probability that they are realized at £ w 3> t* a q is the same as in the true equilibrium. 
In terms of eigen-modes, the above equality is equlivalent to the following one: with 
l* a p defined such that with z(l* a p) ~ (tap)' 1 , 



r a (l) = rp(l) 



for 1 < I < lip. 



(3-4) 
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Fig. 2. The factor fs{l) for different basins 
Fig. 1. The evolution of the indicator R a (tw) (different symbols) around the 6 energy 

for different energy minima a. This is an minima of the system shown in Fig. 1. The 

example of N = 8, T = 0.2 and h = 0.05. error bars represent the scatterings. 



3.2. Slow Modes and Hierarchical Organization of Quasi- Equlibrium Domains 

Let us now introduce what we call slow modes. There are a group of the eigen- 
modes belonging to eigen-values = z(l) < z(2) < . . . z(/ t h) with a threshould value 
lih such that the maxima of their eigen-vectors locate at the elements corresponding 
to the energy minima. We call this set of eigen-modes as slow modes. 



For the slow modes, we have found that equality (3-4) is satisfied within each 
basin B, the neighbourhood of an energy minimum in the phase-space, as 

r«(0 k r B (Z)e M/) for 1< I < l th and a e B (3-5) 

up to some small scattering factors e Sa ^ . In practice, we defined the basin as a 
set of spin-configurations such that T = dynamics starting from any one of it 
converge to the energy minimum with probability one. Then we obtained rjg(l) and 
the scattering e Sa by a x 2 -fitting method. 

The above result means that quasi-equilibrium within each basin is established 
at time scales specified by the slow modes. For instance, we found that the relaxation 
time of the l^-th. mode, which is z^th) -1 , is about 6 MCS at T = 0.2 independent of 
the system sizes we studied (N = 8, . . . , 12). On the other hand, the relaxation time 
of the second mode, which is responsible for the time-reversal symmetry-breaking, 
grows exponentially fast with N, as found previously!*. 

The apparent tree structure recognized in Fig. 1, is due to the follwoing. For 
the factor % of a pair of basins, say B\ and B2, there is a characteristic mode 1% q 2 
such that the equality, 

r Bl il)=r B2 (l) foTl<l<l* BlB2 , (3-6) 

holds within accuracy of the factor e Sa ^ defined above, as shown in Fig. 2. The 
latter means that for t w S> z_1 (^ 1 h 2 )) the two basins are in quasi-equilibrium. Thus 
more and more basins become in quasi-equilibrium with each other as i w increase. 

Let us call the grounp of basins, which are in quasi-equilibrium with each other 
at a given t w , as a quasi-equilibrium domain of age t w . To characterize the growth 
of quasi-equlibrium domain, we measured d max (t w ) which is the maximum of the 
Hamming distance between pairs of energy minima enclosed in a domain of age t w . 
We found that it is an increasing function of i w . 
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solid: response 

broken: from correlation 



3.3. Aging and response to magnetic field 

Lastly, we present a result of our Monte Carlo simulations & performed on 
larger system sizes (but with ,Lj = ±(iV — 1) _1//2 ) simulating the following well 
known experimental procedure HP . For time t w after the temperature quench, the 
system evolves (ages) under zero external magnetic field. Then small magnetic 
field h is swithched on (at t = 0) and the induced magnetization m(t; t w ) is ob- 
served. We also measure the spin auto-correlation under zero-magnetic field, which 
is represented, by means of the notation in the previous sections, as q(t + t w ,t w ) = 

In Fig. 3 we plot the data of m(t; t w ) 
and h[l - q(t + t w ,t w )]/T for N = 
512, T = 0.4 and h = 0.1 against 
t/t w . Interestingly, the two quantities 
coincide with each other, or satisfy the 
flucuation-dissipaton theorm (FDT), in 
the time region t -C t w . The theo- 
rem can only be proven theoretically as- 
suming the true equilibrium situation 
and the above feature seems highly non- 
trival. At present, we speculate that the 
quasi-equilibrium property we discussed 
before is responsible for this: at t <C t w fluctuation and dissipation in the system are 
dominated by those in the quasi-equilibrium domain of age t w , while at later t the 
system starts to invade a bigger domain which causes the breaking of the FDT. 

To summarize, we have studied the relaxational modes and aging of the SK 
model by numerical diagonalization technique and Monte Carlo simulations. We 
focused on the charcteristics of the slow modes and found that the aging proceed by 
hierarchical growth of quasi-equilibrium domains. 

The numerical works have been done on FACOM VPP-500/40 at the Super Com- 
putor Center, Institute for Solid State Physics, University of Tokyo. This work was 
supported by Grand-in-Aid for Scientific Research from the Ministry of Education, 
Science and Culture, Japan. One of the author (HY) was supported by Fellowships 
of the Japan Society for the Promotion of Science for Japanese Junior Scientists. 
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Fig. 3. Response to magnetic field and FDT 
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